Optical analog of Rabi oscillation suppression due to atomic motion 
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The Rabi oscillations of a two-level atom illuminated by a laser on resonance with the atomic 
transition may be suppressed by the atomic motion through averaging or filtering mechanisms. The 
optical analogs of these velocity effects are described. The two atomic levels correspond in the optical 
analogy to orthogonal polarizations of light and the Rabi oscillations to polarization oscillations in 
a medium which is optically active, naturally or due to a magnetic field. In the later case, the two 
orthogonal polarizations could be selected by choosing the orientation of the magnetic field, and 
one of them be filtered out. It is argued that the time-dependent optical polarization oscillations or 
their suppression are observable with current technology. 

PACS numbers: 42.25.Ja, 32.80.-t, 42.50.-p 



I. INTRODUCTION 

The analogies between phenomena occurring in two 
different physical systems open a route to find new ef- 
fects or to translate solution techniques or devices, and 
quite often help to understand both systems better. The 
parallelism between light and atom optics, in particu- 
lar, has been a strong driving force for fundamental and 
applied research, and it has received recently renewed 
impulse with the advent of laser cooling techniques, and 
Bose-Einstein condensation of alkali gases. Atomic inter- 
ferometers or atom lasers are important examples that 
illustrate the fruitfulness of this correspondence. 

The analogy may also enable us to perform in one sys- 
tem experiments which are difficult to carry out in the 
other one. As an example, tunneling time experiments 
are much easier for microwaves than for matter waves 
Q, 0- Typically, optical analogs of quantum mechani- 
cal matter wave effects, if available, are simpler and less 
costly to implement than their matter counterparts. In 
addition, the electromagnetic pulses can be probed in a 
non-invasive way. 

In this paper we shall describe the optical analogs of 
the time-dependent atomic Rabi oscillation in a laser field 
and of several dynamical suppression effects due to quan- 
tum, pure state atomic motion. The Rabi oscillation is 
at the heart of measurement procedures for time and fre- 
quency standards [j| , photon number in a cavity , and 
other metrological applications |5|,|6| , so its dynamic sup- 
pression may be a relevant effect, in particular for atomic 
clocks; This suppression has been also proposed as a way 
to prepare specific internal atomic states by projection 
in quantum information applications 0. 

The relation between light and atom optics is fre- 
quently established at the level of translational (exter- 
nal) degrees of freedom, but here we need, in addition, 
an optical parallel of the two internal levels of the atom, 



which is provided by the orthogonal states of light polar- 
ization. Thus we follow in reverse order, from matter to 
light, a connection that can be traced back historically 
to early experiments of I. Rabi, later modified by N. F. 
Ramsey, which lead to the development of atomic clocks, 
and provided atomic analogs of polarization interferom- 
etry in optics, with the internal states of the atom or 
molecule playing the role of the polarization states of the 
photon |3j. 

For completness a brief review of the dynamical sup- 
pression of the Rabi oscillation is provided in section II, 
which is mostly based on ref. but incorporates also 
some new elements; Section III is devoted to the descrip- 
tion of the optical analog, Section IV presents numerical 
illustrations, and Section V describes a possible imple- 
mentation of polarization filtering making use of mag- 
netic fields. 



II. RABI OSCILLATION SUPPRESSION FOR 
MOVING ATOMS 

Neglecting decay, the effective Hamiltonian for a two 
level atom at rest and in presence of a detuned laser field 
is 
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where f2 is the on- resonance Rabi frequency, 6 = u>l — w 
is the detuning (w^ being the laser frequency and u) the 
atomic frequency) and internal ground and excited states 
are represented as |1) = (*) and |2) = (°), respectively. If 
at time zero the atom is in the ground state, the ground 
and excited components of the state evolved with this 
Hamiltonian arc 
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z l5t/2 [cos{tt't/2) - {iS/Q') sm(fi'i/2)], (2) 

? l5t/2 [{-in/n')sm{n't/2)}, (3) 
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so that the populations ji/n 1 ' 2 ) | 2 alternate oscillating har- 
monically in time with "Rabi period" 2ir/W and effective 
Rabi frequency fi' = (S 2 + iSlj 2 ) 1 / 2 . The oscillation may 
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however be suppressed when the atoms move into a re- 
gion illuminated by a perpendicular laser beam 0. For 
an idealized sharp laser profile in a one dimensional ap- 
proximation (its validity and the three dimensional case 
are examined in 81), the Hamiltonian becomes 



where 



H=p 2 /2m + -Q(z) 
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where z and p are position and momentum operators. 
If reflection is negligible, for moderate to high velocities, 
the stationary wave for incidence in the ground state with 
wavenumber k can be approximated, up to a normaliza- 
tion constant, as 
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This may be interpreted semiclassically as the state of 
an ensemble of independent atoms which travel with mo- 
mentum kh and sustain Rabi oscillations. Note that the 
quantity zm/Tik plays the role of time in the arguments 
of the trigonometric functions, compare with Eqs. @ 
and 0. In other words, the Rabi oscillation is also evi- 
dent spatially in the stationary waves, under the form of 
density undulations of the two components with a "Rabi 
wavelength" \r ; = kh/mVi' . 

The Rabi oscillation may be suppressed "adiabati- 
cally" if a quasi-monochromatic, pure-state wave packet 
enters slowly into the laser illuminated region. Intu- 
itively, and according to a classical picture, the atoms in 
the ensemble start to oscillate at different times because 
of their different entrance instants, so the global oscilla- 
tion averages out if the entrance interval of the packet 
is greater than the Rabi oscillation period. This intu- 
itive interpretation cannot be taken too literally though. 
In particular, the suppression involves a pure state and 
not a statistical mixture. Note also that no incoherent 
"fading" due to decay from the excited state is involved. 
(The effect of fading was discussed in 0.) 

A second group of suppression effects is associated with 
state filtering or velocity splitting. To explain these two 
related concepts, we need a more accurate representation 
than before. In the laser region, let |A+) and |A_) be the 
eigenstates, corresponding to the eigenvalues A±, of the 

matrix \(^» 3>$ ) • O ne easily finds 



A± 
|A±) 



-6±tt' 



1 

-S±Q' b 



(5) 
(6) 



where |A±) have not been normalized. [For later com- 
parison with "circularly polarized states" notice that by 



setting S — and for the case J7=i|fi|, A± 



and 



|A±) 



The stationary state for z > can be 



written as a superposition 

$ fc (z) - C+\\ + )e tk + z + C_|A_}e lfc - z , (7) 
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and the coefficients C+ are obtained from the matching 
conditions at z = The two components (±) have 
two different propagation velocities and in fact the Rabi 
oscillation may be understood as an interference between 
the two terms when k± pa k + m(5 =F O') / \2Tik) . For ini- 
tially quasimonochromatic packets, with the wavenum- 
ber spread much smaller than the average wavenumber, 
Cfc/^o ^ 1; the two different velocities also imply even- 
tually a spatial separation of the |A±) components, so 
that the interference and associated Rabi oscilation fi- 
nally disappear for times larger than the time to split 
the packet into two, 2cr z fc/f2', where a 2 is the spatial vari- 
ance of the wave packet. An extreme case is the complete 
state filtering that occurs for very low kinetic energies, 
when k- is purely imaginary so that e lk ~ z becomes an 
evanescent wave and only the |A+) component survives in 
the laser region for z greater than the penetration length 
[Im(fc_)] . Since the exact form of the surviving |A+) 
state can be modified by the laser detuning and Rabi 
frequency, this state filtering effect provides a projection 
mechanism to prepare especific internal states regardless 
of the incident atomic state Q • 



III. OPTICAL ANALOG 

In the optical analog of the dynamic effects described, 
the internal states will be substituted by orthogonal po- 
larizations of the field, and the laser region by an optically 
active medium; the rotation of the polarization plane and 
corresponding oscillation of the linear polarization inten- 
sities will mimic the Rabi oscillation, and electromagnetic 
pulses will play the role of the atomic wave packets. 1 

The use of a complex electric field for quasi- 
monochromatic pulses within analytic signal theory facil- 
itates the comparison and correspondence between quan- 
tum wave function and electric field components. Some of 
the parallelisms are quite direct. For example, we shall 
keep the same notation for the quantum wavenumbers 
and the optical propagation constants k and k±, or for 
the coefficients in the stationary waves, such as C± and 
the reflection amplitudes R, although their values and de- 
tailed expressions need not be equal. Atomic populations 
in the ground state may be related to total energies in 
the vertical, linear polarization component and, similarly, 



1 Notice also that one further correspondence may be established 
with a spin-polarized electron incident on a region with a per- 
pendicular magnetic field, i.e. there is an analogy between Rabi 
oscillations and Larmor precession, whic h has been used to ex- 
tend the concept of a Larmor clock EjllfJ and other definitions o f 
a traversal time to atomic systems llOt lllll^ [ill ITll 1X^1X11 111 . 
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the atomic excited state will be mimicked by horizontal, 
linear polarization. In our analogy, positions and times 
are present both in optics and quantum mechanics since 
we want to describe and compare pulses and wavepackets 
in space-time. This is at variance with the analogy estab- 
lished by Zaspasskii and Kozlov [18j , in which the spatial 
coordinate played the role of time in a Schrodinger-like 
equation satisfied by the polarization vector. 



A. Basic relations 

From the Maxwell equations in a nonmagnetic, elec- 
trically neutral dielectric medium, the wave equation for 
the E field is 

1 <9 2 E <9 2 P 
Vx ( VxE ) + ^ = -^ ( 9 ) 

where fio is the permeability of the vacuum, c the speed 
of light in vacuum and P the macroscopic polarization 
(volume density of electric dipoles), which we assume to 
be given by the linear constitutive relation 



which, written in components leads to 

Ez = -(X*3 E x + X23 E v)/l3 

and the system 
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E y =0, (17) 



where = 1 + Xjj- It can be solved by making the 
determinant of the coefficients vanish. The result is a 
fourth order equation (k^ + ak^ + (3=0) with coefficients 
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where eo is the permittivity of the vacuum, and x is the 
susceptibility tensor 



X11 X12 X13 
X12 X22 X23 

X*3 X23 X33 



(11) 



The medium is assumed to be non-absorbing so that x 
is hermitian. The absorbing case is briefly considered in 
section Ivl and in the Appendix. 

In the medium, E may not be perpendicular to the 
propagation direction but the electric displacement 



D = e E + P 



(12) 



is perpendicular, since, from Eq. @ and assuming a 
harmonic, plane wave solution e m ' r e~ wt , 



x (k m x E) = -— D. 



(13) 



The subscript m in k m stands for "medium" and it will 
be used to avoid confusion with k = lu/c in vacuum. In 
components, and for k m in z direction, D z = and 
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i.e., the displacement vector components whose ratio de- 
termines the polarization state of the harmonic wave, 
are proportional to the corresponding components of E. 
From Eqs. JTDJl, and (T3), 



k m x (k m x E) 



E 



-XE, 



(15) 



Since we shall assume that the wave incides from a vac- 
uum region adjacent to a semiinfinitc medium, we only 
pick up two physical solutions (if complex they must have 
a positive imaginary part to decay; if real they must be 
positive to implement outgoing boundary conditions) and 
denote them as k± , 



k± = —n±, 
c 



(20) 



n± being the corresponding refraction index. Examples 
will be given soon. 

The polarization corresponding to each solution in the 
x — y plane, up to an intensity constant, may be given 

E x 



in terms of the Jones vector 



E„ 



191. For a closer 



comparison with the atomic case we may use the nota- 
tion E^ = E x , E^ = E y , reminiscent of the atomic 
ground and excited state amplitudes. In particular, the 
internal atomic states ( J ) (ground) and ( ? ) (excited) cor- 



respond to the Jones vectors 



" 1 " 




'0" 







1 



which denote, 



respectively, "vertical" (x-direction) and "horizontal" (y- 
direction) linear polarization. 

Since we are not specifying the total intensity, the po- 
larization is also represented by any proportional vector. 
In particular, it is convenient to work with 



IA+) = 



where 



(Ey/Exh 
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Depending on the value of {^ L )± = a + ib (a and b real), 
the polarization can be linear (b — 0), circular (a = 0, 
b = ±1), or in general elliptic. 

For an optically active medium with a susceptibility 
tensor of the form 



Xn Xi2 
X12 Xn 
X 33 



Re( X 



12; 



0, 



(23) 



the optical propagation constants are given by 

fc± = -(l + Xii±lxi 2 |) 1/2 . (24) 



Assuming Im(xi2) > 0, two orthogonal harmonic solu- 
tions of Eq. J??}, supplemented by Eqs. IjlUfl and for 
right (+) and left (— ) circularly polarized light, forward- 
moving or possibly evanescent (if fc_ becomes imaginary) 
are 



1 
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e + e 
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(25) 



In vacuum, we have instead X = 0, and the following for- 
ward/backward moving orthogonal harmonic solutions 
(with signs +/— respectively) 



(26) 



where k — uj/c > 0. 

If the optically active medium occupies the region 
z > 0, and for normal incidence (k in z direction) of 
vertically polarized light, the harmonic solution reads, 
up to a constant, e~' luJt Fk(z), where 



Ffc(z) =< 



C 4 



(e tkz + R 1 ie- %kz ) + 
' e tk + z + C- 



R 2 ie~ tkz ,z < 



(27) 
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Imposing at z = the continuity of the tangential com- 
ponents of the electric and magnetic fields amounts, for 
normal incidence, to enforce the continuity of the com- 
ponents E XtV and their derivatives, 



1 + iZn 
fc(l-.Rii) 
R21 
-kR 2 i 



Solving the system, 



Rn 
R21 



(fc + fc+)(fc + fc_) : 
—ik(k_ — 

{k + k+){k + k-Y 



k+C+ + fc_C_, 
-iC+ + iC- } 
-ik + C+ + ik-C- 



C+ = 
C- = 
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(28) 
(29) 
(30) 
(31) 



(32) 
(33) 



There are thus two propagation constants k±, and two 
group velocities in the optically active medium, 

du c 
dk± n± 



(34) 



When |xn|i Ixiaj << 1, the two propagation constants 
differ only slightly, 



u 1 + xu±]xi2\ 
c V 2 



(35) 



there is also negligible reflection, Rji » ( j = 1, 2), and 
C± ~ 1/2, so that in the optically active medium, z > 0, 

E (i) M e i(^t) CO s(fc|xi2k/2), (36) 
E (2) w e i(kz- u t) sin (jfc| Xl2 | z / 2 ). (37) 

Each z is thus characterized by some linear polarization 
that rotates when z increases (optical activity). The ro- 
tation length for a full cycle of the moduli squared is 
L = 2irc/uj\xi2\ or, equivalently, the cycle requires a time 
T = 2tt/lo\xi2\ for a quasimonochromatic pulse. Pulses 
are analogous to wave packets in the present correspon- 
dence and are formed by superposition of harmonic com- 
ponents. Moreover, in the quasi-monochromatic regime, 
we can interpret the modulus squared of complex field 
components as (half) short time averages of the real field, 
according to the theory of complex analytic signals [2(J . 
Up to a constant that fixes the actual intensity, the time 
dependent field is given by 



1 



(2tt) 1 /2 



dkA(k)F k (z)e 



(38) 



The time St required, from the entrance instant of the 
pulse peak, to split the pulse into two separated pulses of 
orthogonal circular polarization may be estimated by im- 
posing that the spatial increment between the two com- 
ponents, due to their different group velocities, be equal 
to the pulse width a z . In particular, for xn — and 

|Xl2 1 << 1, 



2at 

|X12 1 



(39) 



where a t = cr z /c. A time dependent observation of the 
polarization rotation, which is the optical analog of a 
temporal Rabi oscillation, requires S t > T to avoid the 
splitting, but also at < T to avoid an averaging sup- 
pression. Combining the two constraints gives the ideal 
conditions 

X12 « ^ « 1. (40) 
k 



IV. NUMERICAL EXAMPLE 

In this section we shall demonstrate with numerical 
examples the time dependence of the polarization rota- 
tion in the optically active medium and several dynam- 
ical suppression effects. The pulse is chosen within the 
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quasi-monochromatic regime, so that we may neglect any 
variation with k of the matrix elements of the suscepti- 
bility and consider constant values in each calculation. 
We have used Eq. (|38l) with the Gaussian amplitude 



A(k) = 



2 

7T 



1/4 



-ikzQ 



(41) 



It is for convenience "normalized" so that J dk \A(k)\ 2 = 
1 and, at time zero, / dz \E^'(z, t — 0)| 2 = 1. In all cases 
the central wavelength is chosen in the visible region of 
the spectrum, Ao = 2n/k = 500 nm. The pulse is thus a 
right moving Gaussian pulse of vertically polarized light 
centered at time t — at zq < 0, outside the optically 
active medium, and with spatial variance cr 2 . 
Figure ^ shows the oscillation of the quantities 



dz\E^(z,t)\ 2 , j = 1,2, 



(42) 



proportional to the total energy in each orthogonal linear 
polarization within the optically active medium versus 
time. For the parameters chosen, the pulse duration is 
much smaller than the rotation period, so the oscillation 
is clearly visible. Whereas optical activity is standardly 
considered in coordinate space and measured in station- 
ary conditions at the end of a slab, here we adopt a dif- 
ferent, time dependent view, closer to the quantum Rabi 
oscillation analog. 




10 20 

t(ps) 



FIG. 1: Time dependent oscillations of the integrated linear 
polarization intensities in the optically active medium. I*- 1 ' 
(solid line) and i' 2 - 1 (dashed line), see Eq. 1421 . are calculated 
for the pulse of Eqs. (1351 and jfij. x = -600 /an, a t = 100 
fs, Ao = 500 nm, X u = 0, xia = 0.0002. 

Complementary perspectives of the time dependent os- 
cillation are provided in Fig. |2 which represent measure- 
ments at fixed positions versus time, or snapshots at fixed 
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FIG. 2: (a) The oscillations of the linear polarization compo- 
nents may also be seen at fixed positions separated by L/2: 
Zl = 1241 /im, z 2 = 2492 fim, z 3 = 3742 fj,m, z 4 = 4993 ;um. 
The pulse arrives at these positions separated by time inter- 
vals T/2, totally dominated, alternatively, by IB^'l 2 (solid 
line) or |i?"'| 2 (dashed line). The parameters are the same 
as in Fig. [I\ (b) Snapshots of the pulse taken at time in- 
tervals T/2: ti = 6200 fs, t 2 = 10368 fs, i 3 = 14536 fs, and 
t4 — 18704 fs. These times correspond to maxima in Fig. Q 
As before, |£ (1) | 2 (solid line) or j£ (2) | 2 (dashed line) dom- 
inate the pulse alternatively at distances separated by L/2. 
Parameters as in Fig. 



times of the intensities of the two orthogonal, linear po- 
larization components in arbitrary units. 

If the pulse width is increased and becomes compara- 
ble or larger than the rotation wavelength, however, the 
oscillation is suppressed, as shown in Fig. [3] This is the 
optical analog of the adiabatic suppression of the Rabi 
oscillation due to a slow entrance of the atoms in the 
laser region. 

Other interesting phenomenon occurs if the observa- 
tion time is large enough so that the pulse splits as a 
consequence of the two different group velocities. This 
is shown first in Fig. 0] where the oscillation eventually 
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FIG. 3: Oscillation suppression due to slow (adiabatic) en- 
trance in the optically active medium. I^ 2 \t) is represented 
for different pulse widths: a t = 100 fs (dotted-dashed line); 
a t = 200 fs (dotted line); a t = 600 fs (solid line). *i 2 = 0.002. 
Other parameters as in Fig. 1. 
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FIG. 5: Modulus squared of polarization field components 
versus time at five different positions: Z\ = 741 /im, 22 = 1341 
fim, Z3 = 1941 fim, Zi = 2541 /im, and z$ = 3141 /im: (a) 
vertical (solid line) and horizontal (dashed line with dots) 
polarizations; (b) right (solid line) and left (dashed line with 
dots) circular polarizations. X12 = 0.08. Other parameters as 
in Fig. 1. 



FIG. 4: Oscillation suppresion of the integrated polarization 
intensities due to pulse splitting: Z^ 1 ' (solid line) and I^ 2 ' 
(dashed line). X12 = 0.08 and other parameters as in Fig. 
1. A transition time that separates the two regimes (with or 
without rotation) may be estimated as the time for the arrival 
of the pulse peak at the medium plus the splitting time of Eq. 
1391 1. This gives 4.5 ps for the present parameters. 

fades away because of the progressive lack of interference 
between the two circularly polarized components. A dif- 
ferent view of this process is given in Fig. [SJi, in which 
vertical and horizontal polarization components are rep- 
resented versus time at five different positions for the 
same pulse of Fig. 0] Notice that each linear polariza- 



tion intensity develops two peaks, to be compared with 
the single peak pulses of Fig. [21 

The gradual separation of the right and left circularly 
polarized components, which for z > are given by 

1 f°° 

£; (±) ee-— / dkA(k)e~ iut C±e ik± *, (43) 
k 1 ' 1 Jo 

is shown explicitly in Fig. 03. The pulse is finally formed 
by two distinct peaks, each of them corresponding to a 
pure right or left circularly polarized subpulsc. 

Finally, the rotation suppression by filtering, analogous 
to atomic state filtering, is illustrated in Fig. El The left 
circular polarization becomes evanescent for the chosen 
susceptibility, xn = an d X12 = 1-2, so, after a tran- 
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FIG. 6: Oscillation suppresion of the integrated polarization 
intensities due to filtering. X12 — 1-2 and other parameters 
as in Fig. 1. In (a) we show the vertical (solid line) and hor- 
izontal (dashed line) polarization components and in (b) the 
right (solid line) and left (dashed line) circular polarization 
components 



sient, the pulse in the active medium is only composed 
by right handed circular polarization, see Fig. |Sp, where 
the total intensities 



dz\E ( 



±)|2 



(44) 



are represented versus time. The absence of left handed 
polarization after the transient peak precludes any os- 
cillation of the linear polarization intensities J^ 1 ' 2 ), as 
shown in Fig. 



V. MAGNETO-OPTIC EFFECTS 

Instead of finding materials with the susceptibility ten- 
sors necessary to observe the different suppression effects, 
it is possible to manipulate \ for an isotropic dielectric 
by applying a static magnetic field B. The susceptibil- 
ity tensor matrix elements in terms of the resonance fre- 
quency wo, plasma frequency u) v , and cyclotron frequen- 
cies uj cu , u = x,y, z, is given in the Appendix for a simple 
Lorentz model. Let us first examine the "Faraday config- 
uration" , with the magnetic field along the z direction. 
In that case the susceptibility takes the form given in Eq. 
(|23|l and the solutions k± correspond to two orthogonal 
circular polarizations. 

If instead the selected magnetic field direccion is x 
(Cotton-Mouton configuration), the form of the suscep- 
tibility tensor becomes 



X = 



Xn o 

X33 X23 
X23 X33 



Re(x 23 ) = 0, 



(45) 



which leads to modes with linear polarizations in x and 
y directions. A field B in an arbitrary direction between 
the Faraday and the Cotton-Mouton configurations pro- 
duces, in general, two elliptic polarizations. Figure [7| 
shows the eccentricities of the two polarizations starting 
from a Faraday configuration and ending up in a Cotton- 
Moutton configuration by increasing the x-component of 
the field. 




FIG. 7: Eccentricity e = \/£ 2 V 2 /^ where £ and 77 are 
the semimajor and semiminor axes of the ellipses associ- 
ated with the two orthogonal polarizations for uj/ljo = 0.98, 

< oj cx /luq < 3.98, uj C y = 0, ui cz /ujo — 2.65. In this case 
(E y /E x )± is purely imaginary, so that the two semiaxes are 

1 and \(E V /E X )±\ 
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A. Evanescence conditions 

Let us determine the conditions which make one, and 
only one, of the propagation constants evanescent, to 
achieve polarization state filtering and selection. We shall 
work out in detail the Faraday configuration, but the 
Cotton-Mouton case or some other magnetic field orien- 
tation could be treated similarly. For the Faraday con- 
figuration, k± are either purely real or purely imaginary, 
and the parameter regions in which one wave becomes 
evanescent are delimited by the zeros of k±, which cor- 
respond, respectively, to 



1 + Xii = TIX12I, 



(46) 



see Eq. (j23J. Using the expressions for xn and X12 given 
in the Appendix, Eq. (|4^|l becomes 



Ul 2 ) 2 — UJ 2 UJ 2 



= ± 



UJ 2 UJUJ C 



[(uj 2 - UJ 2 ) 2 - UJ 2 UJ 2 



,(47) 



with lo cz = eB z /m e . This leads to a second order equa- 
tions for the modulus of uj cz , 

-uj 2 S\uj cz \ 2 ±uj 2 p uj\uj cz \ + S(uj 2 - uj 2 ) \uj 2 - uj 2 ) + UJ 2 p ] = (48) 



where S the sign of [(uj 2 
solutions are 



UJ lc 



4 - H - u 2 ) 



^2cz — 
U4cz — 



(f4-J) 



UJ 

0J 2 p + {uj 2 - UJ 2 ) 



-H-uj 2 ) 



The formal 

(49) 
(50) 
(51) 
(52) 



The parameter region in which filtering occurs purely 
imaginary and k + > 0), is represented in Fig. 8 with a 
light shaded area (in the darker area both solutions are 
evanescent so there is total reflection). It can be divided 
into four w-subregions I, II, III, IV, delimited by the curve 
crossing frequencies lo = luq, uj ~ ujg = (uj 2 , + uj 2 /2) x / 2 , 
and uj = uj 0p = (ujq+uj 2 ) 1 ' 2 . The lower and upper bounds 
for the leftmost region, I, are the curves UJ2 CZ and ujs cz 
(i.e., Eqs. H5Ufl and l|51|) : for region II, the curves uj^ cz and 
^3cz! for region III, uj^ cz and uj^ cz ] and finally the bounds 
for region IV are the curves uj\ cz and uj^ cz . uj = ujs also 
marks a polarization mode change. In the light shaded 
areas, (E y /E x )± = ±i if uj < u>s, whereas {E y /E x )± — 
^fi otherwise. 

An additional constraint is set by the maximum avail- 
able intensity of the magnetic field. It establishes a 
flat upper bound for uj cz , and restricts the frequency 
range where filtering may be accomplished. The smallest 
fields could be used at or near the frequencies ujq and 



2 ; 



3 




0.65 



85 



1.05 

to/co ft 



1.25 



1.45 



FIG. 8: Parameter region in which fc_ is purely imaginary 
and fc+ > (light shaded area) for uj p /loo = 0.84. The criti- 
cal points where curve crossings occur are, from left to right, 
uj/ujq — \ (it separates regions I and II), y/l- + ( lu p/'2cjq) (be- 
tween regions II and III), and \J\ + (uj'p/loq) (between regions 
III and IV). 



In practice the sec- 



would be more useful since uj p could be 



(|49I52|) touch the zero field axis 
ond one, ujq p 

controlled in some cases, e.g. for a gas. Moreover, in 
general the possible perturbation on the filtering effect 
by absorption can be made negligible if the condition 



(uj 2 - 
form 



(uj 2 + uj 2 ) 1 / 2 where the boundary curves in Eqs. 



Jq) » uj~f is satisfied. In ujq p this assumes the 



uj 2 » juj 0p , (53) 



where 7 is the damping constant, see the Appendix. 

Two other important factors that must be taken into 
account for the observability and potential application 
of the polarization filtering are the penetration length of 
the evanescent wave in the optically active medium, and 
the transmittance of the surviving mode. 

The penetration length is proportional to the inverse 
of Im(k_). In practice the medium is not semi-infinite, 
of course, but the filtering effect may be achieved for a 
finite medium as long as it extends beyond the penetra- 
tion length. This should not pose any practical prob- 
lem according to the scales shown in Fig. where 
Ct>o/[clm(£;_)] is shown versus uj cz /ujq for uj = u)q p . 

The "transmission probability" rt+|C+| 2 is the ratio 
between the energy flux of the surviving polarization and 
the inicident flux. It is showed in Fig. ^|for ujq = uj p . 
Up to a significant 25% of transmission may be obtained. 
Finally, real materials hold multiple resonances but our 
analysis can be easily generalized by summing over them 
with appropriate oscillator strenghts factors |19| . 



9 



VI. DISCUSSION 

In summary, we have shown that the atomic time- 
dependent Rabi oscillation and its suppression for mov- 
ing atoms incident on a laser-illuminated region, is anal- 
ogous to the time-dependent polarization rotation and 
its suppression for light pulses entering into an opti- 
cally active medium. The eigenmodes in the laser re- 
gion depend on Rabi frequency and detuning and in 
the optically active medium on the susceptibility matrix 
elements. Since ultrashort femtosecond pulses can be 
tracked experimentally by interferometric photon scan- 
ning tunneling microscopy (PSTM) [2^, the time depen- 
dence of polarization rotation and its suppressions, by 
averaging, pulse splitting or filtering, can be tested ex- 
perimentally. These effects may be of interest for the 
design of an all-optical computer with information en- 
coded, transfered or manipulated using the polarization 
state; similarly, their atomic counterparts may be rele- 
vant in metrology and provide a mechanism for controlled 
quantum internal state preparation by projection or fil- 
tering, i.e., irrespective of the initial state, so that their 
experimental examination at a light-optics level is both 
feasible and worthwhile pursuing. The optical analog of 
atomic state-filtering provides a way to produce polar- 
ized light. Whereas the atomic state may be selected by 
playing with the laser detuning and intensity, the polar- 
ization may be selected by a magnetic field. 




FIG. 10: Transmittance of the surviving mode at to — coop for 

UJp/uJO — 1. 

examine an optical analog of Ref. Q. This work has 
been supported by Ministerio de Ciencia y Tecnologfa 
(BFM2000-0816-C03-03 and HA2002-0002), and UPV- 
EHU (00039.310-13507/2001 and 15968/2004). 
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FIG. 9: Dimensionless penetration length of the evanescent 
wave at uj = u>o p for with lo p /loq = 1. 



APPENDIX A: SUSCEPTIBILITY TENSOR 
ELEMENTS 

The general expressions for the elements of the % ten- 
sor for a homogeneous dielectric in a magnetic field are 
obtained here from a simple, classical, Lorentz model. 
Each electron displacement from its equilibrium position 
r e is assumed to satisfy 



m e 



dt 2 



= -eE - e- 



dVe 

dt 



x B K r e — m e 7 



dr e 

dt 



(Al) 



where K is an elastic force constant that keeps it bound, 
B is a static, external magnetic field, m e the mass of the 
electron and 7 a damping constant. (We neglect the small 
force due to the magnetic field of the optical wave) . As- 
summing that the applied electric field and r e vary har- 
monically as e~ iujt , and using P = —Ner e = xeoE with 
N the number of electrons per unit volume, a lengthy 
but straighforward calculation gives, for 7 = 0, 



X11 
X12 

Xl3 

X22 



(iU 2 — LU 2 ) 2 — UJ 2 UJ : 



V[(OJ 2 -L0 2 ) 2 -U 2 U 2 ](U 2 -U> 2 ) 
2 Lj[itJ cz (ujQ - UJ 2 ) - W cx U> cy w\ 
^[(u; 2 -^) 2 -^ 2 }^ 2 -^) 

2 u)[-iuj cy (u)Q - lu 2 ) - u cx u) cz lv] 



(A2) 



p^I-uj 2 ) 2 -^ 2 ]^ 2 -^) 

I, , 2 , , 2 \ 2 1 i 2 1 > 2 
2 (U) Q —U))—U) LO cy 

^{(ujZ-u; 2 ) 2 -^ 2 }^ 2 -^) 
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X23 
X33 



2 U>[iuJ cx (uj 2 -UJ 2 ) - L) cy U) cz U>} 
0p [(w 2 - W 2 ) 2 - L0 2 L0 2 ](u 2 - U 2 ) 

H( Wo 2 -^)2_^2 ]([J 2_ w2) 



where 



uj = yjK/m e , (A3) 
lo c = eB/m e , lu cu = eB u /m e , u = x,y,z, (A4) 



Ne 2 



m e e 



2 \ 1/2 



(A5) 



are resonance, cyclotron, and "plasma" frequencies, re- 
spectively, and B = (B 2 + B' 2 + B 2 ) 1 / 2 . 



In the hermitian case, i.e., for 7 = 0, Xij 
7^0 the elements in Eq. (IA2D 
form except for the substitution 



would have the same 



K 2 



170;). 



(A6) 



The other non-diagonal elements \ij can be obtained for- 
mally by taking first the complex conjugate of the expres- 
sions of the transpose elements Xji i n Eq. (|A2|) and then 
making the substitution of Eq. I|A6(1 . 
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